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Abstract 

It is an intriguing question how local time can be introduced in the emergent 
picture of spacetime. In this paper, this problem is discussed in the context of 
tensor models. To consistently incorporate local time into tensor models, a rank- 
three tensor model with first class constraints in Hamilton formalism is presented. 
In the limit of usual continuous spaces, the algebra of constraints reproduces that 
of general relativity in Hamilton formalism. While the momentum constraints can 
be realized rather easily by the symmetry of the tensor models, the form of the 
Hamiltonian constraints is strongly limited by the condition of the closure of the 
whole constraint algebra. Thus the Hamiltonian constraints have been determined 
on the assumption that they are local and at most cubic in canonical variables. 
The form of the Hamiltonian constraints has similarity with the Hamiltonian in 
the c < 1 string field theory, but it seems impossible to realize such a constraint 
algebras in the framework of vector or matrix models. Instead these models are 
rather useful as matter theories coupled with the tensor model. In this sense, a 
three-index tensor is the minimum-rank dynamical variable necessary to describe 
gravity in terms of tensor models. 
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1 Introduction 



Tensor models have originally been introduced [IHS] to describe the simplicial quantum gravity 
in general dimensions higher than two, with the hope to extend the success of the matrix 
models in the study of the two-dimensional simplicial quantum gravity. Tensor models have 
later been extended to describe the spin foam and loop quantum gravities by considering Lie- 
group valued indices |1H6]. These models with group indices, called group field theory [7], are 
actively studied with various interesting recent progress [8l-[23]. A key issue is the emergence 
of a new kind of tensor models, called colored tensor models [9], which have more intimate 
topological correspondence to simplicial manifolds than the original versions. There has also 
been a systematic study of tensor models in semi-classical approximations plHST] under the 
interpretation of the rank-three tensor models as theory of dynamical fuzzy spaces [32||33]. and 
the emergence of Euclidean general relativity on emergent spaces has been observed [26l[27] . 

So far the study of tensor models has basically been limited to the cases of Euclidean 
signature. In the standard field theories on fiat spaces, after setting up appropriate causal 
structures, the computations in Minkowski signature can be obtained from those in Euclidean 
signature by means of analytic continuation. However, it is generally not clear how one may 
extend the standard procedures such as Wick rotation to the quantum gravitational situation 
with fluctuating geometries. Moreover, the study of the causal dynamical triangulation [31] 
indicates that the dynamics of quantum gravity in Minkowski signature may be substantially 
different from that in Euclidean signature. The main purpose of this paper is to discuss how 
to incorporate time in the framework of tensor models. 

The advent and correctness of the theory of relativity have established that time is not an 
absolute entity but is rather a relative quantity measured by physical phenomena dubbed as a 
"clock". Since a "clock" is a local object due to the speed limit of light, the definition of time 
is necessarily local and is generally dependent on how the system of "clocks" is organized. The 
principle that physical phenomena themselves should not depend on this kind of ambiguity 
of defining local time provides strong constraints on possible forms of consistent theories of 
nature. 

In the Lagrangian formalism, this ambiguity can well be incorporated by imposing the in- 
variance of theories under the general coordinate transformations of spacetime coordinates. In 
the rank-three tensor models for instance, however, the only dynamical variable is a three- index 
tensor, Mabc {a,b,c = 1, 2, ... , A^), and there are neither space nor locality built in the frame- 
work: a space and its locality are emergent phenomena [35] . Therefore one would have difficul- 
ties in introducing local time into the tensor models and imposing the constraints coming from 
the principle mentioned above. One would try to introduce time t simply as an argument of 
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the tensor as Mabc{t), and impose the invariance of the models under reparametrization t'{t). 
However, this way of introducing time will necessarily allow an entity of a non-local global time 
to exist when a space is emergent. Then emergent field theories on the emergent space will 
not generally satisfy the above principle; emergent field theories on an emergent flat space will 
seriously violate Lorentz symmetry, that is rather strongly constrained experimentally and the- 
oretically [36H38] . Another option to try would be to introduce time for each index of the tensor 
as Mata,,bti„ctc- This option seems to have an appealing feature concerning locality of time. 
However, the contraction of indices of the tensor such as X^afec / dtadtbdtc Mata,bti,ctc Mat.. 
will introduce multiple integrals over times into the Lagrangian formalism. Therefore this op- 
tion seems to require extension of the Lagrangian formalism for multiple time parameters in 
advance, before applying it to the tensor models. 

The above lack of good guiding principles for introducing time to the tensor variable 
would require us to abandon starting with Lagrangian formalism for the purpose. There 
exists another formalism of mechanics, Hamilton formalism, in which time is rather intimately 
related to dynamical evolution of a system. Of course, the two formulations of mechanics are 
equivalent (at least classically), but in the present confusing situation on local time, the latter 
formalism is superior to the former one, because one does not have to know in advance how 
time is represented in the dynamical variable. If necessary, once the Hamilton formalism of 
the tensor models is consistently obtained, one would also be able to obtain the corresponding 
Lagrangian formalism. 

In Hamilton formalism of general relativity, the consistency of dynamics under ambiguous 
choices of local time is guaranteed by a set of first class constraints, which are the generators 
of the local coordinate transformations containing the time direction. In this paper, the set of 
constraints of general relativity will be rewritten in terms of the dynamical variables of tensor 
models, and a rank-three tensor model with first class constraints in Hamilton formalism will 
be presented. This should be equivalent to introducing local time in tensor models. 

The canonical formulations of discrete models of gravity have been discussed in previous 
literatures . An important difference of the present work from the previous approaches 

is that the spatial diffeomorphism is exactly represented by the symmetry of the tensor model. 
Therefore the gauge freedom of spatial diffeomorphism is exactly incorporated by first class 
momentum constraints of the tensor model. Then Hamiltonian constraints will be determined 
by the condition of the closure of the whole set of first class constraints. It is also peculiar 
that time is just a continuous variable unlike some previous approaches [HIHJ]. Thus the 
standard Dirac procedure is applicable, and a rank-three tensor model with local time will be 
formulated in terms of the standard Hamilton formalism with a set of first class constraints. 
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This paper is organized as follows. In Section [21 the general relativity in Hamilton for- 
malism is summarized for the discussions in this paper. In Section [3l the rank-three tensor 
models are briefly overviewed. The limit of usual continous spaces in the tensor models is 
explained. In Section HI an algebra expressed in terms of the canonical variables of a vector 
model is considered, and it is shown that the algebra of the first class constraints of general 
relativity can be reproduced in the limit of usual continuous spaces. In Section [5l however, 
some difficulties in the realization in Section [Hare pointed out, and the necessity of a three- 
index tensor is argued. In Section [6l a rank-three tensor model with first class constraints 
in Hamilton formalism is presented. The momentum constraints represent the kinematical 
symmetry of the rank-three tensor model, and contain the spatial diffeomorphism symmetry 
in the limit of usual continuous spaces. The Hamiltonian constraints are determined by the 
closure of the algebra of the whole first class constraints under the assumption that they are 
at most cubic and respect locality. It turns out to be necessary to break the time-reversal 
symmetry. In Section [TJ it is shown that matter degrees of freedom coupled to the rank-three 
tensor model can be added to the system without destroying the algebraic structure of the first 
class constraints. The matter degrees of freedom can be given by any rank tensors. Section [8] 
is devoted to summary, discussions and future prospects. 



2 The first class constraints from general relativity 



The ADM 



brmulation [H] of general relativity leads to the following algebra of first class 



constraint^ 

{■H{x),n{x')} = e {W{x)6i{x,x') -W{x')6i{x',x)) , (1) 
{Hi{x),n{x')} = n{x)6i{x, x'), (2) 

{ni{x),nj{x')} = ni{x')6j{x,x') + nj{x)6i{x,x'), (3) 

where { , } denotes Poisson bracket, and 6i{x,x') denotes the derivative of the delta function 
with respect to x*. The signature e takes e = 1 and e = — 1 for Minkowski and Euclidean sig- 
natures, respectively. 'H(x) and 'Hi{x) are the super-Hamiltonian and the super- mo mentums, 
respectively, and they are explicitly given by 

4 (.,)^)-^(«<3.-2A). (4) 

Hi{x) = -2Dj7i/, (5) 



*The trivial ones with the conjugate momentums TTNyT^Ni of the lapse and shift variables are omitted. 
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where 7iij{x) is the conjugate momentum to the spatial geometry gij{x), and G and A are the 
gravitational and cosmological constants, respectively. The algebra ([1]), ([2]) and ([3]) can be 
equivalently expressed as 



{H{v), H{w)} = eD{vd'w - wd'v), 
{D{v'),H{w)} = H{v'd^w), 
{D{v'),D{w')} = Div^djw' - w^djv'), 



(6) 
(7) 
(8) 



where H{v) and D{w'^) are defined by 




dx v\x)'Hi{x) , 



dx v{x)'H{x) 



(10) 



(9) 



with v{x) and independent of the canonical variables gij{x) and TTij{x). 

The fundamental roles of the first class constraint algebra ([1]), ([2]) and ([3]) in geometrody- 
namics have been discussed in [17]. An important feature is that the right-hand side of ([1]) 
contains the inverse spatial metric g'^^{x) to raise the index of 'Hi{x) to W{x). Therefore the 
constraint algebra has structure functions depending on the canonical variable, and is not a 
Lie algebra with constant structure constants. To reconcile the constraint algebra with the 
spacetime diffeomorphism, the on-shell conditions "H = "Hj = must be imposed. In addition, 
under some (reasonable) assumptions, the form of the super-Hamiltonian has been shown to 
be given uniquely by the form (jl]). It is also argued that the constraint algebra does not 
change by adding matters with non-derivative couplings with gravity. 

3 A brief overview of tensor models 

This section will give a brief overview of the rank-three tensor models (with no time) to prepare 
for the discussions in the following sections. As explained in Section [T], tensor models have 
some variations with distinct interpretations. This paper deals with the rank-three tensor 
models, which have a three-index tensor as their only dynamical variable. The rank-three 
tensor models can be regarded as theory of dynamical fuzzy spaces [^IB3] . This interpretation 
of the tensor models is especially convenient to understand the semi-classical behavior of the 
tensor models: classical solutions are regarded as background fuzzy spaces which approximate 
continuum spaces, and the perturbations around the classical solutions are interpreted by 
effective field theories on the spaces. In fact, various classical solutions and the perturbations 
around them are studied to show the phenomena of emergent Euclidean general relativity 
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on emergent spaces [23H^ . It is especially noteworthy that, as theory of dynamical fuzzy 
spaces, the rank-three tensor models can deal with any dimensional spaces, unlike the original 
proposals [TW3], in which the ranks of tensors are related with dimensions. 

Let me denote the three-index dynamical tensor by Mabc- The tensor is assumed to satisfy 
the generalized Hermiticity condition [IHS], 

M,,, = Mb,, = M,,, = M,;, = M:,, = M;,,, (11) 

where * denotes complex conjugation, and the indices run as a,b,c = 1,2, . . . , N. Because 
of the generalized Hermiticity condition (fTTj) . the symmetry which can be associated to the 
rank-three tensor models is the orthogonal group 0{N), 

M:,, = 0/0,''0/Ma>t'c', Oa' E 0(iV), (12) 

instead of a unitary group of an hermitian matrix model. 

While a continuous manifold can be described by a coordinate system, a fuzzy space is 
defined by an algebra of functions on it. The algebra of functions (pa {a = 1,2, . . . , N) can be 
characterized by its structure constants Cabc aiO 

4>a4>b = Cabc4>c- (13) 

While a usual continuous space can be characterized by a commutative and associative algebra 
of functions, a noncommutative space for instance can be characterized by a noncommutative 
associative algebra. One may even consider a nonassociative algebra to define a nonassociative 
space [HHSI]. 

To relate fuzzy spaces to the configurations of the rank-three tensor models, let me intro- 
duce an inner product [32] . 

(0a|06) = Sab, (14) 

which is assumed to be bi-linear. Now let me identify the structure constants with the dy- 
namical variable of the tensor model, 

Cabc = Mabc. (15) 

This physically means that the rank-three tensor models are interpreted as theory describing 
the dynamics of fuzzy spaces. The identification (ITSl) and the generalized Hermiticity condition 
(|TT|) imply the following cyclicity property on the algebraic structure, 

{(Pa(pb\4>c) = {4>a\4>b4>c) = {4>b\4>c4>a) , (16) 



Repeated indices are assumed to be summed over throughout this paper, unless otherwise stated. 
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and also the properties on complex conjugation, 

0a = (17) 
(0a06)* = Ma. (18) 

These properties ([TB]) . ffTTl) and (ITS]) characterize the fuzzy spaces which can be associated with 
the configurations of the rank-three tensor models. The fuzzy spaces with these properties have 
various interesting properties concerning symmetries, uncertainties, and reduction procedures 

[321IS2H51]. 

In the analysis of emergent Euclidean general relativity from the rank-three tensor models 
[21H2n], the following particular form of Mabc with Gaussian functions has extremely been 
useful. The indices of functions (pa are assumed to be given by the coordinates of a usual 
continuous D-dimensional space, 

a = x = (x^a;^...,x^), e R, (19) 

and Mx^x2x:i is assumed to be given by the following Gaussian form, 

M^.,., = Bg{xif'^g{x2fl^g{x^f'^e^^ X2)' + rf(x2, X3)' + rf(x3, Xi)^)] , (20) 

where B and /3 are positive constants. The fuzziness of the spaces is characterized by the 
parameter /3. Here a metric tensor field gij{x) is assumed to exist on the D- dimensional 
space, g{x) = Det[gij{x)], and d{xi,X2) denotes the geodesic distance between two points Xi 
and X2- The form of (1201) respects the diffeomorphism symmetry, as g{xY^^ guarantees the 
diffeomorphism invariance of an index contraction: M^^^^^M^^^^^ = J dx^,/g{x) ■ ■ ■ . 

The Gaussian configuration (|20l) is merely an idealized working hypothesis which singles 
out the modes corresponding to those of general relativity. This hypothesis has very well 
explained the qualitative features of the results of the numerical analysis [251429] . which have 
shown the emergence of Euclidean general relativity on emergent spaces. The detailed values 
of Mabc inay be different from that, because the actual degrees of freedom of the rank-three 
tensor models are discrete and finite, while they are continuous and infinite for the continuum 
index (1191) . Also the simple Gaussian damping form is idealizing M^-^^j-^rc,^ which is locally 
distributed with respect to the relative locations of Xi. Thus the Gaussian configuration fl20l) 
should be regarded as an infrared effective idealized description, which would be obtained 
from a coarse-graining procedure [21]. 

The algebra ([13]) with ([15]) and ([20]), 

0^10^2 = ^Sx2X3</'x3, (21) 
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defines a fuzzy space of the kind satisfying f fTB]) . f lT7|) and ffTS]) . Intuitively, the function (^j. 
represents a fuzzy point at x, which has fuzziness with a length scale ~ Fo^' the fiat 

case gij{x) = 6ij, the Gaussian configuration fl20l) becomes 

Mitx3 = ^exp ((xi - xa)^ + {x^ - x,Y + {x^ - x,f)] . (22) 

This configuration indeed respects the Poincare symmetry and can be considered to represent 
a flat D-dimensional fuzzy space [oT| . In the limit /3 — )■ oo with an appropriate normalization 
B, the fuzzy space algebra approaches 

(pxi(px2 = S^^i^i - X2) (pxi in /3 00. (23) 

This is the limit to the usual continuous spaces with no fuzziness. Indeed the algebra (!23l) 
is commutative and associative. In the rest of this paper, this limit will be denoted by the 
pointwise limit, and will be used to reproduce the constraint algebra of general relativity in 
Hamilton formalism from tensor models. 

In fact, in the subsequent discussions, the Gaussian form fl^Ul) or fl2^ is not essential. It is 
merely a representative of the configurations of Mabc, which have relatively local distributions 
when the 0{N) gauge symmetry f|T2|) is appropriately fixed, and have a parameter that can 
be tuned to take the pointwise limit f l2^ . Another important thing is that one cannot take 
the pointwise limit as a starting point. As will be seen, one has to start with a finite fuzziness 
and then take the pointwise limit. For simplicity, I will only use the fiat expression (122 p in 
the subsequent computations, but it should be straightforward to extend to the general cases 
by using the diffeomorphism invariant expression (!20|) . 



4 Realization of constraints by a vector model 

In this section, I will try to incorporate the constraint algebra of general relativity ([1]), ([2]) and 
(|3]) in the framework of vector models in Hamilton formalism. The discussions will proceed 
almost well, but some difficulties, which will be discussed in Section will arise in a final 
step. Matrix models will also be abandoned due to the same difficulties. 

The degrees of freedom of a vector model in Hamilton formalism are assumed to be given 
by Ma (a = 1, 2, ... , A^) and their conjugate momentums tTq. They are assumed to satisfy the 
canonical relation of Poisson bracket, 

{Ma,7T,} = 6ab. (24) 
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Let me define 



(25) 



■V 



(26) 



[a6] 



where e' = ±1 is a signature, and its relation with e in ([T]) or ([6]) will be given later. The 
C^'s satisfy C^^^^ = Cj^^^ and C^^^ = —C^a]^ respectively. As will be explained in detail be- 
low, (125 p and fl2Bl) mimic the super-Hamiltonian and super-momentums of general relativity, 
respectively. 

From f l24p . f l25p and (1251) . one can straighforwardly obtain the following algebraic relations. 



which look similar to (P), (IZ]) and of general relativity. Here Hy and Dy are defined by 



and the upper indices S and A indicate the symmetric properties of the matrices, f ff, = v^^ and 
v^fj = —v^jj, respectively. The v,w's are assumed to be independent of the canonical variables. 
The square bracket [v, w] denotes the commutator of matrices, [v, w]ab = VacWcb — WacVcb- Thus 
C^^j^j and form a Lie algebra under the Poisson bracket. Especially, due to ( l29i) and the 
anti-symmetry of f^^,, form the Lie algebra so{N) of the orthogonal group, which is the 
kinematical symmetry of the vector model. 

In the following, I will show how the constraint algebra of general relativity (E]), (JTj) and 
([8]) can be reproduced from the algebra of a vector model (l27j) . (I28l) and (l29l) . Let me assume 
that a situation similar to the Gaussian configurations in Section [3] is occurring in the vector 
model. Then the indices are assumed to take the coordinates of a continuous Z)-dimensional 
space as in (|T9|) . 

Let me first discuss (12^ . In (12^ . the computation of the Poisson bracket has been reduced 
to a commutator of matrices. Let me consider a (infinite-dimensional) matrix in the form. 



{Hv{v''),Hv{w'')} 

{Dy{v^),Hy{w')} 
{Dy{v^),Dy{w^)} 



-e'Dyi[v',w']) 

Hy{[v^,w']), 
Dy{[v'',W^]), 



(27) 
(28) 
(29) 



Hy{v')^v!,Cl,^, 

Dy{v^)^V^,Cl^,^, 



(30) 
(31) 



v: 



xy 



(32) 
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which is ant i- symmetric v^y = —Vy^. Here f*(a;) are assumed to be arbitrary smooth functions 
on the Z)- dimensional space. Then the commutator between two such matrices is given by 

[v'^,w\y = dz {v\x) + v\z))5^ix, z){w\z) + {y))5j{z,y) - (x ^ y). (33) 

Because of the derivatives of the delta functions, the computation of the right-hand side of 
(I33p tends to become cumbersome. It is much more straightforward and easier to do the 
computation by considering a test function f{y). By multiplying the right-hand side of fl33|) 
with f{y) and performing the partial integrations over y and z, one obtains 

dzdyfiy) [{v\x) + v'{z))5i{x, z){w^ (z) + {y))Sj{z,y) - (x ^ y)] 
= 2v'{x){didjW^{x))f{x) + 4v'{x){diW^{x))djf{x) - {v ^ w). (34) 

Comparing the last expression of fl5^ with the right-hand side of 

^ dy{v\x)+v\y))6,{x,y)f{y) = {d,v\x))f{x) + 2v\x)dj{x), (35) 
one obtains 

[v^,w\y = ^ {[v,w]\x) + [v,wny)) 6,{x,y), (36) 

where 

[v,w]\x) = v^{x)djw\x) - w^{x)d,jv\x). (37) 

Comparing fl5Bl) and flTTI) with fl52]) . one concludes that fl2^ with fl5^ exactly reproduces 
([8]). Thus Dv{v^) is the analogue to the generators of spatial diffeomorphism in the general 
relativity. 

In the next place, let me discuss (1251) . Consider 



< = c(/3) w y-^) exp(-/3(x - yY), (38) 

where {x—yY = {x—yy{x—y)\ The ( !38|) satisfies the symmetry w^y = Wy^. The function w{x) 
is assumed to be an arbitrary smooth function on the D-dimensional space and its argument 
in (l38l) takes the middle point between x and y. The coefficient c(/3) depending on (3 will be 
determined later. 

The Gaussian form in (138|) follows the Gaussian configurations in Section [3l and, at the 
final step of the following computations, the pointwise limit /3 — )■ oo will be taken to compare 
with the constraint algebra of general relativity. This Gaussian form is considered just because 
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of its simplicity. In fact, the following discussions do not depend on the details of the form. 
What is needed is that has distributions within finite ranges of relative distances between 
X and ?/, and one can finally take a smooth pointwise limit. 

The computation of the Poisson bracket (!28|) has been reduced to the commutator of the 
matrices ( 15^ and which is given by 



xy 



""^^^ ^ dz {y\x) + v\z))5i{x, z) w I ) exp(-/3(^ - y f) + {x ^ y) 



c(/3) 



div\x)w ( 1^ ^ ) exp(— — y)^) 



+ 2v'{x)d^ [ w [ ) exp(-/3(x - y)') 



(x^y). (39) 



The last expression is rather confusing, because the derivative with respect to x in the last line 
produces a factor /3, which makes the pointwise limit difficult to handle. It is again much easier 
to do the computations by considering a test function f{y). By multiplying the right-hand 
side of fl5Ul) with f{y) and integrating over one obtains, in the leading order of 



dyf{y) 



c(/3) 



div\x)w [ ^^-^ ] exp(-/3(x - yf) 



+ 2v\x)d^ ( w ( ) exp(-/3(a; - y)') ] + {x ^ y) 



c{(5) / rf;zexp(-/3z2) [t;*(x)9,w(x) + 0(/3-i)] f{x). 



(40) 



This concludes 



X + y 



exp (-/3(x - yf) 



(41) 



in the leading order of 1//3. Thus, by comparing the right-hand side of (l4T]) with ( l38l) . (l28l) 
with f p2|) and fl38|) exactly reproduces ([7]) in the pointwise hmit (3 oo. 

Finally let me discuss ( !27l) . What should be computed is the commutator between the 
matrices in the form f p8|) : 



v^,w%y = cilSy I dz 



X + z\ f z + y 
V I w 



{x ^ y) 



exp [-p{^x-zf-f3{z-yf). 

(42) 



To systematically do the computation, let me again consider a test function /(y) and evaluate 



c{(5f / dydzfiy) 



X + z\ ( z + y 
V I — - — w 



{x ^ y) 



exp(-/3(a;-^)2-/3(z-2/)2). (43) 
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By a change of variables, y = a + b + x, z = b + x, the integration becomes 

c{(3f j dadbf{a + b + x) [v{x + b/2)w{a/2 + b + x)-{v ^w)] exp(-/3a^ - pb'^). (44) 

Because of the Gaussian damping factor for a and 6, the large /3 limit of this integral can 
well be evaluated after Taylor expanding in a and b the integrand other than the exponential. 
Then one obtains 

\2n-D-l I 1 / J^^2 ^2■^ 



c(/3)'/3-^-^ (^-^ J dzz' exp{-z') j dzexp{-z' 

dif{x) {v{x)diw{x) — w{x)div{x)) + ij^f{^) {v{x)didiw{x) — w{x)didiv{x)) . (45) 
Comparing with ( 135|1 . this concludes 

[v',w%y = cic(/3)2/3-^-i^ {k{x) + k{y)) 5,{x,y), (46) 
where ci is a numerical factor, and 

li{x) = v{x)diw{x) — w{x)div{x). (47) 



The ugly contraction of the indices in fl46|) comes from the simplified assumption z = z^z^ 
in fl55]) . If a general metric is assumed as z"^ = gijZ^z\ one will have 

j da^ a'a^ exp(-/3a2) ~ I3'^~^g'^, (48) 

and P6l) will contain . It would be straightforward to do the computations in a full covariant 
fashion following the expressions in Section [3l 

Thus, by taking 

c(/3) = /3'^q^ (49) 

e' = -e, (50) 



one obtains ([6]) of general relativity from ([2j 

An important fact in the above computations is that the finiteness of the range \x — y\'^ < 
of the distribution in f l38|) plays an essential role in the derivation of fj^6|) . even though the 
pointwise limit /3 — )■ oo is finally taken. This is indicated by the extra factor in ( l46l) . This 
is also the reason why the final expression depends on the inverse metric g^^ , which appears 
also in the constraint algebra (EI) of general relativity. If the matrix (l38ll were assumed to have 
a full diagonal expression like 5{x,y) for instance, ([6]) would not be reproduced. 
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5 Difficulties and a solution 



From the discussions in Section HI the super-Hamiltonians in the vector model should be given 
by C^ab) ( |25|) . Then the generators of an infinitesimal local "time" translation will be given 
by Stiab)C}^ab)- Howcver, in the regime discussed in Section [31 the indices label each "point" in 
a space, and the "time" is generally highly non-local, since x and y in can freely take 

any values. Therefore this "time" is very different from the conventional notion of time in 
physics. One may instead consider a "time" in a diagonal form 6tab ~ ^ia^ab- But, as discussed 
in the last paragraph in Section HI such a full diagonal form cannot correctly reproduce the 
constraint algebra of general relativity. 

There exists also another more general and serious difficulty. Under the Poisson bracket, 
the constraints fl25p and f l26p form a Lie algebra with constant structure constants. This means 
that the future time evolution is completely determined by the action of the Lie group element 
parameterized by the "time" . As reviewed in Section [21 this is substantially different from the 
gravitational case, in which the constraint algebra has structure functions depending on the 
metric. One cannot expect non-trivial dynamics to occur from the constraint algebra of the 
vector model. 

These problems cannot be solved by considering a matrix model. One would try a set of 
constraints, 

Cfab) = T^acT^bc + e'M,,Mfe„ (51) 

Cjl] = 7r,,Mb, - M,,7rb,. (52) 

But the problems above appear in the same as in the vector model. 

A solution can be given, if there exists a three-index tensor Mabc- By using Mabci the two 
indices of C,\ab) "^^^ contracted as 

Cl^M,,,Cl,y (53) 

Then the infinitesimal time can have only one index as btgC^. Moreover, in the regime 
explained in Section [3l SixM^yz will provide a non-diagonal distribution of a finite range for 
y, z, that is necessary to reproduce the constraint algebra of general relativity in the pointwise 
limit /3 — )■ oo, as is discussed in the last paragraph of Section [H And also, as will be seen in 
Section [21 the constraint algebra does not have structure constants, but they rather depend 
on Mabc and its conjugate. 

In principle, the indices of the super-momentums of the vector model CJ^^ can also be 
contracted by Mabc- However, this is not a valid option, because, as explained in Section [H 
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the tensor models have the orthogonal group symmetry, and Cj^^^j are the generators of this 
kinematical symmetry. If they were contracted, the number of constraints would be reduced, 
and the symmetry could not be fully incorporated by the constraints. 

Another important thing is that Mabc must be a dynamical variable. If not, the Poisson 
bracket between the super-momentums Cj^^,] and would not close: 

{viCl.^^v^M^.^Cj,,^} = viv^M,,,{Cl,^XU} (54) 

= [v^,vMUCl,^, (55) 

where vM denotes a matrix vMab = v'^Mcab- The expression in the last line does not in 
general have the form of linear combinations of f lS^ . Therefore one has to include Mate as a 
canonical variable, and make it transform appropriately under the Poisson bracket with the 
super-momentums. 

The discussions in this section can be summarized as follows. To correctly reproduce the 
constraint algebra of general relativity and the usual notion of time, it is necessary to include a 
three- index tensor Mabc as a dynamical variable. In this sense, Mabc is the dynamical variable 
corresponding to gravity. The vector and matrix variables may be incorporated consistently 
in the constraint algebra, but they are not necessary. Instead they can rather be regarded 
as some matter degrees of freedom, which can be added consistently, as will be discussed in 
Section [71 



6 Realization by a rank-three tensor model 

The discussions in Section pimply that the pure gravitational system should be obtained from 
a rank-three tensor model which has a three-index tensor as its only dynamical variable. The 
canonical variables are assumed to be given by Mabc and iTabc- They are assumed to satisfy 
the generalized Hermiticity condition fITT]) and its conjugate correspondence, 

T^abc = T^bca = T^cab = T^bac — '^acb — '^cba^ (56) 

respectively. The Poisson bracket between the canonical variables is assumed to be given by 

{Mabc, T^dcf} = Sabcdef = ^ad^be^cf + ^achf^cd + ^af^bd^cc- (57) 



Let me consider 



C(ab) = 2 {T^acdT^bdc " ^MacdMbdc) , (58) 
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Clab] = - {T^acdMbcd " MacdT^bcd) , (59) 

where the relation e' = — e in f lSOp has aheady been used. Let me define 

H'{v'')^v!,C^,b), (60) 
Div^)^vtAat], (61) 

where ffj, = and f^f, = — f^, and they are assumed to be independent of the canonical 
variables. Then H' and D satisfy 

{H'{v^), H'{w^)} = eD{[v^, w^]), (62) 
{D{v^),H\w')} = H\[v^,w']), (63) 
{D(^^),D(^^)} = D([t;^^^]), (64) 

which are actually the same as ( 1271) . ( 128|) and ( 129|) of the vector model. 

The discussions in Section H] for a vector model uses only the fact that the algebra of 
constraints under Poisson bracket can be written as commutators of matrices as in f l27p . 
and fl29|) . Since this is the same for the tensor model as in fl62|) . f l63|) and f lM|) . the constraint 
algebra of general relativity can be reproduced in the same way as in the vector model. 

As discussed in Section [5l let me contract the indices of C(ab) with Mabc to construct the 
possible super-Hamiltonians defined by 

Ca = Mabc{C{bc) + A(5{,c), (65) 

where I have included a new term with a coefficient A. This term is meaningful because it 
adds a non-constant term to Ca- The term is consistent with the kinematical symmetry of the 
tensor models, since 6ab is an invariant of 0{N). 

Now let me define 

H{V) = VaCa, (66) 

for an infinitesimal parameter Va independent of the canonical variables. Because of the 0{N) 
invariant form of f l66|) . it is obvious that 

{D{v^),H{w)} = H{v'^w), (67) 

where v^Wa = v^j^Wb, and therefore the Poisson bracket between D and H closes. 
The Poisson bracket between two H^s is given by 

{H{v), H{w)} = {Cab, Ccd}VabWcd + {Vab, Ccd}WcdCab + {Cab, Wcd}VabCcd + {Vab, Wcd}CabCcd, 

(68) 
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where I have defined 



Cab = C(^ab) + A(5ab, 
Vab = VcM(^ab)c, 
Wab = WcM(^ab)c, 



(69) 
(70) 
(71) 



for notational simphcity. Here M(^ab)c (and 'iT(^ab)c) is the symmetrization defined by 



M(^ab)c = - {Make + M^ac) , 



(72) 



7r(afe)c = - (TTafec + T^bac) ■ 



(73) 



Note that, thanks to ( ITTi) and ( l56|) . M(^ab)c and 7r((jfe)c are totally symmetric with respect to 
their indices. The computation of the first term in (1681) is similar to ( 162|) . The last term in 
(|68|) trivially vanishes. The sum of the second and the third terms can be computed as 



The first term of the right-hand side is what is desired. The matrices v and w are distributed 
by Mabc as in flTOj) and flTTl) . Therefore, this term exactly reproduces the constraint algebra of 
general relativity in the pointwise limit /3 — )■ oo, as has been discussed in Section HI However, 
the second term is proportional to neither C[ab] nor Ca- Therefore Ca are not appropriate as 
super-Hamiltonians, since the constraint algebra does not close. 

To cancel the second term in (1751) . let me first proceed under the assumption that the 
super-Hamiltonians be invariant under the time-reversal transformation T^abc —T^abc- I also 
assume that the terms are at most cubic and local: each term must be connected, not like 
MabbMcdeMcde ^OT instaucc. Then the possible cubic terms other than those in (^5^ can be 
listed as 



Ve{Mi^ab)e, T^cfgT^dgf - eMcfgMdgf}WhM^cd)hCab - {v ^ w) 
= Ve7^cfgS{ab)e,dgfWhM(^cd)hCab - {v ^ w) 
= 2 (VaWb - WaVb) C^d 'r^{ac)e M(M)e- 



(74) 



Thus the result of ( 168|1 is obtained as 



{H{v), H{w)} = eD{[v, w]) + 2{VaWb - WaVb) Ccd Vr(ac)e M(^bd)e- 



(75) 



= T^abcT^bdeMa 
Cf^ = Mabcn 
Cf^ = TTabcMbcdT^d 



cdei 



lee) 



(76) 
(77) 
(78) 
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(79) 

(80) 



where I have not cared the order of the indices, since this is not essential in the following 
discussions. It can easily be shown that C^'^\ . . . , C^^^ are not useful to cancel the second term 
in (1751) . When C^^^ is added to (1651) . the Poisson bracket will have additional contributions, 
{C^^\C^^^} and {C^^\C}. From exphcit computations, {C^^\C^^^} is proportional to super- 
momentums, and might be allowed. However, {C^^^ , C} contains terms 

VaWb{C^^\Cb} ~ di VaWeMdefTlacdMfijMcji + (^2 VaWbMdefMbgeMgcfKacd H , (81) 

where dj's are some integers. The first term in (1811) can be used to cancel part of the unwanted 
term in (1751) . but the second term in flHTj) is problematic. This term is not proportional to 
the super-momentums, and it can also be shown that this term cannot be canceled by adding 
C'^'^\ . . . ,C^^\ Thus it is not possible to construct super-Hamiltonians that contain (^Si), are 
cubic at highest and local, and respect the time-reversal symmetry. 

A simple solution can be found, if one relaxes the time-reversal symmetry, while the other 
assumptions are kept intact. Let me consider 

Ca = (Mabc + e'nabc) (C^bc) + XSbc) , (82) 

where e" is a parameter which will be determined in the following. This will become the final 
form of super-Hamiltonians. Define 

Hiv) = VaCa (83) 

for an infinitesimal parameter Va independent of the canonical variables. Because of the 0{N) 
invariant form of (1831) . one obtains 

{D{v^),H{w)} = H{v^w), (84) 

where v'^Wa = v^b'^b- 

The Poisson bracket between iJ's is given by 

{H{v), H{w)} = {Cab, Ccd}v'abW'^d + {^afc. Ced}</ab + {Cab, w'^dV^b^cd + {^fe, w'^^}CabCcd, 

(85) 

where 

v'ab = 'Vc{.M^^ab)c + (^"T^(ab)c), (86) 
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(87) 



Computations similar to the case of (!68|) result in 

{H{V),H{W)} = eD{[v',w'])+2{l-e{ey) {VaWb - WaVb) C,d7l(ac)e M^M)e. (88) 

For the closure of the constraint algebra, e" is determined as 

e" = ±l/^e. (89) 

Then 

{H{v),H{w)} = eD{[v',w']). (90) 

The commutator in fl90|) is of the matrices fl86l) and flHTj) . This contains not only Mq^c 
but also Tiabc- But, as is stressed in Section HI the constraint algebra of general relativity can 
be reproduced without the details of the distribution of these matrices. One would be able 
to expect that, in the regime discussed in Section [31 not only Mabc but also the conjugate 
momenta TCabc have relatively local distributions which become pointwise in the limit /3 — t- oo. 
Therefore the addition of iiabc will not affect the derivation of the constraint algebra of general 
relativity in the pointwise limit. 

In the Euclidean case e = — 1, fl89l) implies that e" must take an imaginary value e" = ±i. 
Then an infinitesimal time-like shift will generate 

SMabc = {Mabc, H{V)} = ±tVaCbc + " • • . (91) 

This direct appearance of the imaginary unit i will generally violate the generalized Hermiticity 
condition (|TT|) . Therefore the Euclidean case is not consistent. 



7 Coupling with matters 

Although the vector model considered in Section H] is not suited for describing gravity, it may 
be added as a matter. Because of the kinematical character, the total super-momentums are 
simply given by the sum of those of the three-tensor (159!) and the vector (l26l) as 

C total — /7 I /^V fc\0\ 
[ab] = <-lab] + <^[ab], [^^) 

which are obviously the generators of the 0{N) transformation. The total super-Hamiltonians 
would be defined by a summation, 

^totai ^^^^ (^^^^ ^ e"7r,,,)C(l) , (93) 
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where Ca and C^^^^ are the super-Hamiltonians of the three-tensor (1521) and the vector fl2^ . 
respectively. Here the coupling between gravity and matter is realized by the contraction of 
^(ab) with Mabc + ^"T^abc- To chcck whether this actually produces the desired result, let me 
compute the Poisson bracket, 

(94) 

where 

H'°"'\v) = VaC'^K (95) 

The first and the last terms in flM|) produce the desired form, and the fourth term obviously 
vanishes. The computation of the sum of the second and the third terms is essentially the 
same as (185|) . and one obtains 

if the condition (|89|) is satisfied, where 

^*°*'^'(^'^) = ^a1,C&f, (97) 

and v' and are defined in (!86|) and (!87|) . respectively. Because of the apparent O(A^) 
invariant forms, it is straightforward to derive 

It is obvious that one can add various matters in similar manners. It will be straightforward 
to extend fl2^ . fl25l) . fl26|) . fl92l) and fl93|l to such general cases. Any rank of tensors and 
statistics, bosonic or fermionic, will be allowed. 



8 Summary, discussions and future prospects 

In this paper, I have discussed how local time can be introduced in tensor models. Since it 
was not clear how to formulate such tensor models in the Lagrangian formalism, a rank-three 
tensor model with first class constraints in Hamilton formalism has been presented. The 
discussions have made it clear that a three-index tensor has a prominent feature necessary 
for the purpose. The other rank tensors can be added as matter sectors which has coupling 
with the three-index tensor. In this sense, the rank-three tensor model, which contains a 
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three-index tensor as its only dynamical variable, can be regarded as the gravitational sector 
which can universally couple with matters. 

The momentum constraints have straightforwardly been constructed so that they incorpo- 
rate the kinematical symmetry of the rank-three tensor models. Then the consistency of the 
local time evolution requires Hamiltonian constraints and the momentum constraints to com- 
pose a closed first class constraint algebra. This closure condition gives strong limitations on 
the possible forms of super-Hamiltonians, and they have been determined on the assumptions 
that they are local and cubic at highest. They have turned out to contain terms which break 
the time-reversal symmetry. It has also been shown that the constraint algebra reproduces the 
first class constraint algebra of general relativity in the limit of usual continuous spaces. The 
first class constraint algebra closes exactly without any approximations, and this is a good 
feature which marks distinction from the previous approaches to discrete gravity in 

the canonical formalism. 

The most important implication of this paper is that it is in principle possible for the 
rank-three tensor model to have the first class constraint algebra which reproduces that of 
general relativity in the limit of usual continuous spaces. However, the limiting procedure is 
rather formal assuming the regime discussed in Section [31 Therefore it will be necessary to 
investigate whether the limit can be generated by the actual dynamics of the rank-three tensor 
model. Not only the constraint algebra, but it would also be necessary to check whether the 
equations of motion of general relativity can directly be reproduced in the continuum limit 
from the super-Hamiltonians. Especially, the fate of the terms with the breaking of the time- 
reversal symmetry has to be studied. If the breaking effect has turned out to be too large to 
be allowed by experiments/observations, the super-Hamiltonians must be reconsidered under 
relaxed assumptions. By canonical transformations of the variables, it is also possible to 
change the property of the super-Hamiltonians under the time-reversal transformation. 

Regardless of whether such a regime in Section 12] actually exists in the dynamics, the 
canonical rank-three tensor model of this paper has a significance of its own as a model of 
emergent space and gravity with consistently incorporated local time. The form of the super- 
Hamiltonians has similarity with the Hamiltonian of the c < 1 string field theory , which 
has applications to random surfaces. Since the c < 1 string field theory may be regarded as 
(l-l-l)-dimensional gravity with matters, the rank-three tensor model of this paper may have 
some connections to it. The connections may provide an interesting testing ground as well as 
some hints for the dynamics of the tensor model. 
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